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g^ : ABSTRACT 

By analysing the infinite-dimensional midisuperspace of spherically sym- 
D , metric dust matter universes, and applying it to collapsing dust stars, one 

finds that the general quantum state is a bound state. This leads to a dis- 
crete spectrum. To an outside observer, the geometry is static if the initial 
radius of the collapsing star is smaller then the Schwarzschild radius. In that 
case the discrete spectrum implies Bekenstein area quantization: the area 
of the black hole is an integer multiple of the Planck area. Knowing the 
microscopic (quantum) states, we suggest a microscopic interpretation of the 
thermodynamics of black holes: by calculating the degeneracy of the quan- 
tum states forming a black hole, one gets the Bekenstein-Hawking entropy 
(the entropy is proportional to the surface area of the black hole). All other 
thermodynamical quantities can be derived by using the standard definitions. 
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1 Introduction 



Recently a simple minisupcispace describing the Oppenheimer - Snyder (OS) 
collapsing star was found [1]. The semiclassical wave function of that model 
(e.g. the semiclassical sohition of the Whcclcr-DcWitt equation) is a bound 
state. This leads to quantization conditions. The corresponding (Bohr- 
Sommerfeld) quantization condition can be written in the form 

F(M,i?o) = ^(n+^) , n = 0,l,2,... (1.1) 

where M is the mass of the collapsing star, Rq is its initial radius, and 
F(M, Ro) is a function of M and Ro to be given later [1]. For fixed i?o, (1-1) 
implies mass quantization. 

The idea of mass quantization is an old one. Using general arguments 
from quantum mechanics (of adiabatic variables) and general relativity, Beken- 
stein [2] got the black hole area quantization condition 

Ml^^hn , 71=1,2,3,... (1.2) 

where Mj^ is the irreducible mass of the black hole, which is related to its 
surface area , A — lOvrM^y . This discrete spectrum can be related to the 
thermodynamic properties of black holes [3]. 

Wc will show that in the case of black holes, one can get (1.2) from 

(1.1) , which is correct for any OS star. Beside the quantization condition 

(1.2) , one can find in this explicit model also the general quantum states 
(the solutions of the Wheeler-DeWitt equation) forming this star. By that 
one can hope to understand better not just the quantization conditions, but 
also the thermodynamical properties of black holes. This is the purpose of 
this paper. 
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While it has been known for two decades, there are still some open ques- 
tions concerning the entropy of black holes. The classical considerations [4], 
which use the analogy between some geometrical properties of black holes and 
thermodynamics, gives the "laws of black hole thermodynamics" , from which 
one can get the entropy. Semiclassical considerations [5] , on the other hand, 
use the (formal) path integrals (approach) to find the partition function, and 
then the entropy. But both approaches do not use basic statistical mechanics 
reasoning, namely, finding the entropy by calculating the number of differ- 
ent "microscopic states" that correspond to the same "macroscopic state" 
that we call a black hole. This "missing hnk" is very important, because 
it requires the understanding of microscopic states forming a (macroscopic) 
black hole. Those microscopic states are the quantum states, so their impor- 
tance to the understanding of black holes thermodynamics is obvious. The 
OS model gives a one dimensional minisuperspace, in which of course one 
cannot hope to get a degeneracy that will give a nonzero entropy. So one 
must extend the model. One such an extension is the inclusion of inhomo- 
geneous (spherically symmetric dust) distributions. This was done a long 
time ago by Lund [6]. Lund used the dust matter as a "clock" and then 
fixed the gauge completely, reducing the constrained matter-gravity theory 
to an unconstrained one. We will use Lund's infinite dimensional "midisuper- 
space" , and apply it to the collapsing star case. Though infinite dimensional, 
Lund's midisuperspace shares some resemblance with the OS model, so one 
can analyse it in the same manner and find the quantum states that corre- 
spond to a classical black hole. Knowing the microscopic (quantum) states, 
one can calculate their degeneracy, and find the entropy of the black hole. 
Then by using the standard thermodynamical definitions one can get the 
other thermodynamical quantities (e.g. the temperature). 
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We will consider both the static "eternal black hole picture" [7], and the 
dynamic Hawking evaporation one [8]. They both can be studied in our 
framework, and the results that we get are in agreement with the known 
ones. 

In this work we use the semiclassical approximation only. This is for 
two reasons: first, the OS model as well as Lund's one, are correct only 
semiclassicaly. And second, we use Einstein gravity (coupled to matter) 
which should be (at least) a good approximation semiclassicaly. 

We use geometrical units G — c — 1. 

The paper is organized as follows: in chapter 2 we describe the OS model, 
solve (semiclassicaly) the corresponding Wheller-DeWitt equation, and get 
the mass and area quantization conditions. In chapter 3 we describe Lund's 
midisuperspace and find the (semiclassical) solutions to the Wheeler-DeWitt 
equation. In chapter 4 we describe the midisuperspace of a collapsing star, 
and find the quantum states forming this star (black hole). In chapter 5 we 
study the thermodynamical properties of black holes in this framework, and 
chapter 6 presents some concluding remarks. 



2 The OS Model 

2.1 The OS Minisuperspace 

In 1939 Oppenheimer and Snyder [9] found a very simple solution (of Einstein 
gravity couple to dust matter) describing a collapse of a spherically symmet- 
ric homogeneous dust star. In their solution the Schwarzschild exterior is 
smoothly connected to the interior region, which is a slice of a Friedmann 
universe. 



4 



The interior region is described by the Friedmann hne element 



ds^ = -N'^{t)df + a^{t)[dx^ + sin\dnl 



(2.1) 



The range of x is < x < xo , where Xo < 7r/2 . At x = Xo the interior is 
matched to the exterior Schwarzschild solution. If M and Rq are the mass 
and initial radius of the star, then the matching conditions are 



M = -aosin^xo 
Ro = oosmxo 



(2.2) 



where flo is the initial Friedmann radius. 

The gravitational Lagrangian may be split into its interior and exterior 
parts, 



fXO 

Lq^A-k sin xdx 
Jo 



^d'-3Na 



+ 



f 

Jr>rs 



gRd'^x 



(2.3) 



where is the surface radius of the collapsing star. 
The matter Lagrangian is 



X 



(2.4) 



where p is the density of the star and is the four-velocity of the matter 
particles. Energy momentum conservation, V^T^*^ = 0, implies p = po/a^, 
where po is a constant to be determined by the initial conditions. The OS 
model requires po — Sao/Sn. So using (2. 3), (2. 4) and C/''C/^ = — 1 we get the 
total Lagrangian 

^2 



Lg + L 



M 



rxo 
Jo 



sin^X 



aa 



— N(a — ao) 



+ [ ^/^Rd^x (2.5) 

Jr>rs 



The Hamiltonian corresponding to (2.5) is 

1 



Aa^a 



P„ + Q;o(a - oo) 



+ / Hd'^x 

'r>ra 



(2.6) 



where Pa = dL/dd, and ao = 127r J^'^ sin^xdx ■ Because the classical solution 
for r > is the Schwarzschild space-time, for which Rsch. = 0, only the first 
term in (2.5) (or (2.6)) will contribute to the semiclassical dynamics^]. 

The Wheeler-DeWitt equation is the quantum version of the classical 
Hamiltonian constraint, OH /ON = , in the coordinate representation: 
|\E' >= ip{a) , Pa = —ihd/da . Using (2.6) we get the "Schrodinger equation" 

("^'^ + V{a)j Ha) = (2.7) 

where V{a) = Aala^a — ao) . If we define x = a — ao/2, we get 

I -^'^ + j ij{x) = Ei){x) (2.8) 

where uj = 4ao and E = agOQ. As we can see from (2.8), ip describes an 
harmonic oscillator (with mass m = 1/2). So the semiclassical wave function 
describes a one-dimensional harmonic oscillator, which is of course a bound 
state. 

As in the Hartle-Hawking case [10], the solution of (2.8) describes a su- 
perposition of two "universes''^, one that collapsses to form a black hole, and 
one that expands, a white hole. 



^Using the path integral approach, the semiclassical wave function is 

"^wKB = Aexp[iSciass./TT] 

and we see that only the first term in (2.5) will contribute. 
^ The solution can be written as 

where p{x) = y/\V{x')\dx' . 
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2.2 Mass and Area Quantization 

Because the wave function of the OS model describes a bound state, the 
spectrum is quantized. Semiclassically we should use the Bohr-Sommerfeld 
quantization condition, but in the case of an harmonic oscillator it is exact, 

E{n) = nuj(n + -] , ra = 0,1,2,... (2.9) 



2 

Using the definition of E and uj (see below (2.8)) we have 

^aoal = h(n+^^ (2.10) 

In [1] we consider only the case Rq >> 2M which corresponds to the usual 
cosmological situation. In that case we have ao — 47rxQ, and using (2.2) and 
(2.10) we get [1] 

MRlin) = ^h'(n + ^^ (2.11) 

For fixed initial radius, (2.11) gives mass quantization. 

Of course all the above describes a dynamical process: the collapse of the 
star. In this work we try to understand the quantum properties of a "static 
black hole," as seen by an outside observer, and especially to find its entropy 
and temperature (as measured by that observer). For an outside observer, 
the above picture can be static only if Rq < 2M . In that case the geometry 
outside the horizon is always Schwarzschild, which is of course static. The 
OS model requires Rq > 2M [9], so only if Rq = 2M this model can describe 
a static geometry everywhere outside the horizon. 

In the case Rq = 2M we have Xo = ^/2, and we get from (2.2) and (2.10) 

M\n) = ^h(n + ^^ , (2.12) 

a surface area quantization: the surface area of the black holeQ goes linearly 
with the quantum number n. 

^ For a Schwarzschild black hole Mir ^ M,&o WttAP. 
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This is Bekenstein's result [2], but we got it by using an explicit model, 
and by solving the corresponding Wheeler-DeWitt equation. The fact that 
we use a very simple (and even a non- realistic) model (the OS star), and still 
get his results, suggests that these are quite general. 

In Bekenstein's original paper [2] the prefactor for was h/2 , but 
further considerations by Mukhanov [3] suggest that the prefactor should 
be hln2/4n. In our case the prefactor is h/STf^. This prefactor is a model 
dependent, and the best that we can hope (using our simple modelQ) is to 
get the same order of magnitude. This is in fact what we got. 

In the OS model, the (microscopic) state of a black hole with a (macro- 
scopic) mass M, is |\E'„ >, where n satisfy (2.12). So in this model, for each 
macroscopic state (labeled by the mass M, or the area A) there is only one 
microscopic stateQ (labeled by the quantum number n). So the entropy of the 
black hole in the OS model is zeroQ. This is because the OS minisuperspace is 
"too small" (one dimensional). If we want to understand the thermodynam- 
ics of black holes, we must extend the model. In the next chapter we describe 
a much bigger midisuperspace (an infinite dimensional one), which will turn 
out to be much more appropriate for studying black holes thermodynamics. 



3 Lund's Midisuperspace 



^ For example, we take only Rq — 2M. A more reasonable model should take some 
average between Rq — Q and i?o = 2Af. 

^ We have a one dimensional harmonic oscillator, so there is no degeneracy. 
''The entropy goes like 

S ~ /n(number of microscopic states) = ln{l) = 
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3.1 The ADM Reduction Process 



For a global hyperbolic space-time, M — R x T,^^\ one can use the ADM 
splitting [11], and write the Hamiltonian of a dust matter coupled to Einstein 
gravity in the form [6] 



H 



I 



drxdt 



N{n'^ + £) + Ni{W + V') 



(3.1) 



where and Aj are the lapse function and shift vector respectively. and 
W are the gravitational super-hamiltonian and super-momentum 

1 



«J 2 



TT 



(3.2) 
(3.3) 



where hij is the induced 3-metric on T.^^^ [h = det{hij) ), and tt*-' its conju- 
gate momenta (vr = 7r*j ). S and are the dust matter Hamiltonian and 
momentum respectively 



(3.4) 
(3.5) 



where is a unit normal vector to the hypersurface E^^) , and T^^, is (using 
(2.4)) 

T^, = 167rpC/^C/, (3.6) 
If we define the scalar field 6 , 



(3.7) 



we can treat it as a dynamical variable (and p, which is not a dynamical 
variable, will be a function of hij, tt'-', 0, V^^^, to be determined later.). 
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In the spherically symmetric case one can use the {R, 9, 0)-coordinates on 
E*^^^ in which 

dsl^) = e^^dR^ + e^^dnl (3.8) 

where /i and A are functions of t and R, and dQ| is the volume element in 
5"^. If we define tt^ and tta such that 

•1 o„ 1 ox 1 



n'^ = diag (^^e-'^TT^, ^e-^V,, -e-'^'sin-'eTrx^ (3.9) 



we get 



g-2M |^2A" - 2AV + 3(A')')]) (3.10) 

^ _e-2/^(7r;-/iV^-AVA) (3.11) 
£ = (l67rp/iV2)-^^2 (3 12) 

V"" = p^/i^V (3.13) 

where prime denote differentiation with respect to R, and = dL/d(f) — 
— IdnNh^/'^U^p. We see that if we choose the coordinates for which — 0, 
and using U^U^ = -1, we get p = (167rh^/^y\l + h^^{(t)'f)-^/'^p^ , so from 
(3.12) we get 

£^{l + h'\<P')'fp^ (3.14) 

and wc sec that £ goes linearly with p^. This suggests that we can use as a 
time variable. And indeed taking (p = —t , = gives (using (3.10)-(3.13)) 
the known general solutions [12]. 

To complete the gauge fixing (the reduction process) one must choose also 
the i?-coordinate. From the equations of motion (derive from (3.10)-(3.13)) 
one can get that X'e'^~^ is a function of R only, so one can choose 

R = aV"/^ (3.15) 
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Now the conjugate momenta are 

TTR = -{X'e^-^y\, (3.16) 
TT, = 7r,-e^((A')-^e-V)' (3-17) 

Using (3.11) and (3.13), the supermomentum constraint, H^ + V^ — 0, is 
now 

TTR + TTxX' +P4>(j)' = (3.18) 

After solving the constraints, one ends up with the reduced Lagrangian (or 
Hamiltonian) 

Sred ^^'^ j dtdR{'Kyy - Hadm) (3.19) 



where 



and 



y = 8e^ (3.20) 
TT, = ^e-^ (3.21) 



Hadm = i?' f Jtt^ + {2R)-\R-^ - l)y\ (3.22) 



So we end up with oo^ unconstraint degrees of freedom, the y{r) field, with 
the Hamiltonian (3.22). The space of all 7/(r)-field solutions is what we call 
"Lund's midisuperspace" . 



3.2 Quantum States 

The reduced Hamiltonian is 



H = 47: J HADMdR = 4tt J R^ (^^tt^ + {2R)-\R-^ - l)y^ dR (3. 



23) 
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We use the coordinate representation (in Lund's midisuperspace) 



y = y (3.24) 

h 5 , ^ 

TXy = -— 3.25 

I oy 

So the corresponding Schrodinger equation (in a different context we could 
call it the Wheeler-DeWitt equation, see sec. 2) is 

«^ . / , 1 _ (3.26) 

where |\E' > is the quantum state, and \E'[?/;t] =< ^\y{R,t) > is the wave 
functional. In this representation it is a functional of the field y{r) and a 
function of time t. As one can see, we use the "yvrj^-ordering"! in (3.26), but 
different ordering will not change our results, which anyway are correct only 
semiclassicaly. 

A very important feature of (3.23) is that there are no /^-derivatives in 
H. This meens that the infinite degrees of freedom (d.o.f.) are decoupled. 
Let Rs be the surface of the dust ball, then we can divide Rs to equal 
parts, Rk = {k = 1,2, N). The "continuum limit" is N oo. The 
vector space {y{r)} is now {y = (?/i, ?/2, 1/Af)} where yk = y{Rk), and the 
Schrodinger equation (3.26) becomes 



d<f{y,t) 



N 



9t h 



' Vk dyl 



^{y,t) (3.27) 



where = Rl and bk = (-R^^ l)/4 are positive constants {bk is positive 
because < i? < 1 [6]). Because of the decoupling we can write |\E^ > as a 
direct product 

1^ >= > 1^2 > ■ ■ -l^jv > (3.28) 
®The field y is always to the left of its conjugate momenta. 
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and from (3.27) we have now 



ih- 



dt 



-ctk—TT^ + hUk 



Vkdyl 



(3.29) 



where ^kiVk^t) —< ^k\yk{t) >■ The ADM Hamiltonian is time-independent 
so we have Ti-ADMiVk) — Ek — const, and 



-ak 



Vk dy] 



2 + hVk 



If we define 



<lfk{yk:t) ^ Ek-^kiVk) 



1 



Xk = yk- -Ek 



we get the following harmonic oscillator Schrodinger equation 

+ n^kl^l^l "^kiXk) = ^k'^kiXk) 



(3.30) 



(3.31) 



(3.32) 



y 2mk dx\ ' 2 

where ruk — l/2ak , ouk — V^dkH a-nd = mk{ukE}.)^ The solutions 
of (3.32) are >= \nk > , where jn^ > is a one-dimensional harmonic 
oscillator with energy = hcuk{nk -|- 1/2) , or 

^al^'El ^h(nk + I) (3.33) 

So the space of quantum states describing this spherically symmetric dust 
"universe" is spanned by 

{\^ni,n2,...,nN >= 1^1 > \n2 > • ' >JV >} (3.34) 



The total energy is 

N N 



E = Y.Ek = Y. [miRk' - Rf)Vf\nk + (3.35) 

k=l k=l 

In the N ^ oo limit, one can have a finite energy only if one uses the Wick 
order, so -|- 1/2 must be replaced with n^, but we will come to that later. 
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4 Midisuperspace for a Collapsing Star 
4.1 The Homogeneous Case 

In the homogeneous case (the OS case), it is convenient to use (2.1). The 
3-metric in the {R, 6, 0)-coordinates, is 

/ijj = diag ^a^ ^^^^ , a^sir?Xi siri^xsir}^^^ (4.1) 

and from (3. 8), (3. 15) and (4.1) we get 

R = X'e^-P = cosx (4.2) 

We see that R is not a "usual radial coordinate". For example the origin 
(r = 0) is = 1, and Rmin = cosxo > 0. From (3.15) one can see that 
R = Tc/ve [6], where is the circumference radius, and r^. is the extrinsic 
radius of curvature. Because the horizon is a minimal area, R{horizon) = 0. 
One can get this explicitly in the OS model, because a static black hole 
corresponds to xo = ^/2, so R{horizon) = cosxo = 0. For asymptotically 
flat space R goes to unity at spatial infinity. It is convenient to use a different 
radial coordinate]^ 

r = sinx = (1 - R^f^^ (4.3) 
Now one should replace (3.23) with 

/32rVl - r2p2 

67ir 



H 



r=„,/N, 3/:rV i - -T" 37rr \ , ,,,, 

/ n{r)dr= [—^ -+ o n 2 ^r^ 

^0 JO \ 37r y 1\J\ — I 



where y = 8e^ = 8ar, and Py = "iTiyy/QAr^/l — r^. It is very easy to see 
that (4.4) is the correct Hamiltonian, because from (2.6) and (4.4) we have 
®The coordinate r grows with the usual radial coordinate while R decreases. 
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Py = 37rrP„/2ao\/l 



— r^, so 



sinxo 



= 127r / dr , — ^ + a 

p2 

— + aoa (4.5) 



4Q;oa 

which is exactly the gravitational part of (2.6). 
The Schrodinger equation is 

^1* >= > (4.6) 

In Lund's formalism E is an undefined constant, but in the collapsing star 
case, we must impose the matching conditions: the solution must be smoothly 
matched to the outside Schwarzschild space-time. This will constrain E. 
From (2.1) and (2.4) we have 

E = 327T^ / sin^xdxa^P (4-7) 
Jo 

Because E is time independent, we can easily calculate it at the beginning 
of the collapse, when the star is at rest. Using a{t = 0) ~ Qq and p{t = 0) = 
3M/AnRl, we get 

E^?^ (4.8) 
Now we can use the matching conditions (2.2) to get 

E = aoao (4.9) 

and the Schrodinger equation (4.6) is exactly (2. 7), and one can get the results 
of chapter 2. 

Notice that we could get (4.9) from the requirement that the collapse 
start {t = 0) at rest. We will use that in the inhomogeneous case. 



15 



In the homogeneous case the difference between Lund's midisuperspace 
and a collapsing star midisuperspace is that in the latter the energy E is 
constraint. In the inhomogeneous case the constraints are more complicated, 
but one can deal with them in a similar way. 



4.2 The Inhomogeneous Case 

We saw that (in the homogeneous case) one can use Lund's formalism, impose 
the energy condition (4.9), and get the OS results. This can be generalized 
to the inhomogeneous case. 

Let us first go beck to the homogeneous case. Using (4.4) and (4.9), the 
Schrodinger equation (4.6) can be written as 



32rVl - r2 Snr 



37r y 2vl - r 



^[y]c/r = (4.10) 



where yo = y{r,t = 0) = 8aQr. 

Now we can generalized to the inhomogeneous case. The result (3.23) (or 
(4.4)) are correct for any spherically symmetric dust ball (not just for the 
homogeneous one). In particular, they are correct in the case of a collapsing 
(spherically symmetric dust) star. If we want the collapse to start (t = 0) at 
rest, which is the generalization of (4.9), then we must haveQ 

S{r) = :^^^=yo (4.11) 



2VT 



r 



because only in that case y{r,t = 0) = 0. So the form of (4.10) is quite 
general: it is correct also in the inhomogeneous case. The only difference is 
that in the inhomogeneous case, the Friedmann radius, a, can be a function 
of r too. So as a function of (r, t) the field solution is different, but it 
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E = £{r)dr, so n{r) = £{r), and from (4.4) and (4.10) we get (4.11) 
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has the same form y{r,t) = 8a{r,t)r. The space of all field solutions that 
describe a collapsing star is of course a subspace of all the field solutions 
(Lund's midisuperspace) , but at this stage we do not need to know the specific 
restrictions; what is important is that we can use (4.10). 

After making the discretization rk — rgk/N {k — 1,2, N) , we get from 
(4.10) 

E + Pkivk - yf ) m = (4.12) 

k=l \ Vk ) 

where = (32rjky'l - r^)/37r , = Zirrk/I^l - rl and yf^ = yo{rk). Using 
^{y) — Ilk ^k{yk) we have a set of N independent equations 

[at^ + (3k{yk - yf )) *fe(yfe) = (4.13) 

Defining Xk = yk — yk^ /'^ , we get the (harmonic oscillator) Schrodinger 
equation 

+ -rrikcolxl ^^kixk) = ek'^kixk) (4.14) 



V 2mkdkl 2 

where = l/2ak , ujk — Sr^ and ek — uikCoKy^^) /8. The quantization 
conditions are 



.(0).2 . / 1 



rkrukiyD =hi^nk + -) (4.15) 
The total energy is 

N (0) 

^-37rE (4.16) 



k=i 2\/l-rl 



20^ 

and from (4.15) we get 



N 

E^hJ2nk{nk + l/2) (4.17) 

k=l 



where Q,k — 4:/aj^\ i^k^ — ^('^fe, t — 0)) 
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So a quantum state describing a collapsing dust star, starting at rest, can 
be written as 

l^(star) >= |ni > |n2 > • • -Iun > (4-18) 

where jn^ > is a one dimensional harmonic oscillator (exited to the level Uk) 
with frequency cuk — 8rfe. Each \nk > and so |^ > are bound states, and we 
end up with quantization conditions. 

In the homogeneous case y^^ — SaoTk , (a^f^ = oq) and from (4.17) we get 

4?? 

£^w = — EK + V2) (4.19) 

and using (4.9) we have 

laoal = hY,{nk + l/2) (4.20) 

k 

which is (2.10) (remember that in the homogeneous case there is only one 
d.o.f. so N = k = l) . 

In the general inhomogeneous case, fi^ is not /c-independent, but we can 
use "mean field" reasoning and write (4.17) as 

N N 

E = nJ2 ^fcK + i/2) = h<n>J2{nk + 1/2) (4.21) 

k=l k=l 

One can use (4.21) as a definition of < Q >, which is the "average" of Qk- 

The results (4.18) and (4.21) are correct for any collapsing spheriacly 
symmetric dust star, which start at rest. But in the case of a static black 
hole, one can relate E and < Q > to the mass of the black hole. In the 
homogeneous case we have Cq = 2Mf,h so from (4.9) we have E = Ott^M^^, 
and from (4.17) < Q >= Qk = (2Mft/j)~^. In the inhomogeneous case this 
can be generalized by dimensional arguments to 

E i Mbh (4.22) 
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wher ~ denote equality up to a constant of order unity. Now using (4.21) 
and (4.22) we get 

TV 

Mlh -^^T.i^k + 1/2) = h{ntot + 1/2) (4.23) 
fc=i 

Which is the generaUzation of (2.12) to the general inhomogeneous case, in 
agreement with Bekenstein's result. 

5 Black Hole Thermodynamics 
5.1 Entropy 

Using the standard statistical mechanics definition, the entropy of a (macro- 
scopic) black hole is 

S,H = ln[A4,(M)] (5.1) 

where we choose to be one, and J^bh{M) is the number of microscopic 
states that correspond to the same macroscopic state (a black hole) with 
mass M. The microscopic states are (4.18), and for a given M (or a given 
E), J\fbh{M) is the number of different |*(star) >-states, that satisfy (4.21) 
(or (4.23)). 

Consider first the hmit N ^ oo: according to (4.21) the energy, E, will 
be finite only if we use the Wick order. Then if only a finite number of d.o.f. 
are excited, E will be finite. In that case there are infinitely many other d.o.f. 
that are in their ground states, and we face two (probably related) problems: 
first, our semiclassical approximation is very bad when most of the d.o.f. are 
in their ground state. Second, N'bh is infinite (there are infinitely many ways 
to choose a finite number of exited d.o.f. from an infinite number of total 
d.o.f.) so the entropy will diverge. As a matter of fact, the limit — > oo is 



19 



questionable. For example (2.10) is as much a radius quantization condition 
as it is a mass quantization condition. This means that one cannot simply 
divide infinitely many times. This is clearly a quantum gravitational issue. 
But there is a way to avoid it: If we choose not to use the Wick order, then 
if E is finite, must be finite too. We see that (2.12) or its generalization 
(4.23) provide us with a natural cutoff^ Nmax ~ {M/MpY, which for a 
classical black hole is a big, but still finite. 

Of course 1 < < Nmax- In a sense, N describe the amount of in- 
homogenity. For = 1 we have the homogeneous case, for N = 2 the 
"almost" homogeneous one, and so on until = N^ax which describe the 
general inhomogeneous star. We have then 

^ max 

Ar,h{M) = J2 ^bh{N,M) (5.2) 

where Mbh{N,M) is the corresponding number for a specific A^. Our semi- 
classical approximation is good as long as A^ is much smaller then Nmax, but 
we will see that the contribution to (5.2) from A^ > Nmaxl'^ is the same as 
from A^ < Nmaxl'^- So at least we have a good estimate to (5.2). 

It is easy to see that the A^'s that will contribute to (5.2) must satisfy 

N = Nmax - 2j , j = 0, 1, {Nmax " l)/2. Let US start from A^ = Nmax- 

In that case we have only one state (4.18), >= |0 >i |0 >2 ...|0 >N,nax 
. Next we consider A^ = Nmax — 2, it is easy to see that they are Nmax — 2 
statesQ |\I' >= |0 >i ..|1 >k --\0 >Ar„^^-2- In a similar way, we have for any 

A^ = Nmax - 2j 

AUiV, M) = (5.3) 

^^In our geometrical units Mp = Ip = h^^^ . 

^^The oscilators are distinguishable because they have different frequencies, ujk = Sr^. 
Or in other words: different k's correspond to different shells which are distinguishable 
because they have different radii. 
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where are the binomial coefficients. So 

N,n{M)= Y: (5-4) 

j=0 

Though there is no known analytic expresion for (5.4) [13], it is elementary 
to show numericaly that 

n 

Yl ^L-j - exp(0.962n - 0.320) (5.5) 

j=0 

In our case we have n = {N^ax — l)/2 ~ {M/MpY, so 

A4/.(M) ~ exp (c^^ (5.6) 

where C is a constant of order unity. And we find (using (5.1)) that the 
entropy of the black hole is 

S,H -Cj^ + So (5.7) 

The entropy is proportional to the surface area of the black hole, or equiva- 
lently, it goes linearly with the quantum number ritot, see (4.23), in agreement 
with the Bekenstein-Hawking entropy. 

We could use a different approach to calculate the entropy: One can use 
the Wick order, so (4.22) is replaced with 

N 

E^h<n>J2nk (5.8) 

k=l 

but still take a finite Nmax- From (2.10) and Rq ~ oq we have that ARq > 
h/Ro, and from Rq — a^rs we get (Ait!o)mm — O'oT^s/Nmax- So in the case 
of black holes, = 1 and i?o ~ M , we have N^^ax ~ Rl/^ ~ M^/M^. 
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In this case N can take all the integer values between N^ax and unity; the 
degeneracy is 

KniM)^ E = 2^--' (5-9) 

7V=1 

Now the entropy is 

Sbh^ln2^ + So (5.10) 

We see that (5.7) and (5.10) have the same form, the entropy is proportional 
to the surface area of the black hole, but the pref actors are different. One 
should get the correct prefactor in the full exact model. 

Another thing to notice is that the degeneracy (and so the entropy) of 
the gravitational d.o.f. is very similar to the degeneracy of other field d.o.f. 
[14], so it is tempting to think that the gravitational d.o.f. that we use are 
more appropriate for a unified scheme. 



5.2 Temperature 

Using the standard thermodynamical definition of the temperature we have 

= (5.n, 

We have Ef,fi ~ M, and using (5.6) we get 

(5.12) 

in agreement with the Hawking temperature [8]. 

One may argue that the microscopic states (4.18) are microscopic both 
to a freely faUing observer ( "Kruskal observer" ) and to an outside observer 
( "Schwarzschild observer"). This means that (5.1) should be the same for 
Kruskal observers as well. Then the entropy (5.7), and temperature (5.12), 
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are the same for both the Kruskal and Schwarzschild observers. This con- 
tradicts the known results that a thermal Schwarzschild state corresponds 
to a zero Kruskal temperature [7,14]. But remember that though the micro- 
scopic states (4.18) are the same for both observers, the macroscopic states 
are quite different. For a Schwarzschild observer there is an horizon, and 
from the no-hair theorems, there is only one macroscopic quantity (in the 
Schwarzschild case) by which one can determine the state. This is of course 
the mass M of the star. In that case the degeneracy is exactly what we get 
from (4.21), and indeed we have (5.7) and (5.12). On the other hand, for 
a Kruskal observer, there is no horizon, and the macroscopic state is deter- 
mine by an infinite number of macroscopic quantities. For example in our 
model, we have a global hyperbolic space-time, so a classical solution is de- 
termined by the initial data0 {y{r,t = 0),y(r, t = 0)). In our case we have 
y{r,t = 0) = 0, so a classical solution is determine by y{r,t = 0), or equiv- 
alently by all the moments P„ = Jq" r"'y{r,t = 0)dr, which are macroscopic 
quantities. This means that all the fi^'s in (4.23) are determined, and (at 
least semiclassically) the state (4.18) is determined completely hj the macro- 
scopic state. This means that for a Kruskal observer there is no degeneracy, 
and the entropy and temperature vanish. 

The entropy (5.1) is sometimes called "entangeled entropy", but we think 
that (at least in the case of black holes) it should be consistent with the 
thermal entropy. The way to check this is to couple the system to other 
fields. When we couple the gravitational d.o.f. to other fields, we have 
the following picture: The fields are in "equilibrium" with the gravitational 
d.o.f. (the black hole). According to a Schwarzschild observer, it is a thermal 
equilibrium with the temperature (5.12), and according to a Kruskal observer 

^•^The question of classical (and of course quantum) observables in gravity is an open 
one, but in principle one should be able to determine those quantities. 
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it is a zero temperature situation. This is a static scenario, in agreement with 
the "eternal black hole" picture [7], and with the Thermo-field approach [15]. 

One can consider also a dynamical situation: a black hole creation and 
evaporation. This will be done in the next subsection. 

5.3 Hawking Evaporation 

So far we have studied only a static picture as seen by a Schwarzschild ob- 
server, in which a black hole is in a thermal state in equilibrium with the 
outside region. But one can use our formalism to study also the dynamical 
process of Hawking evaporation. Using the semiclassical adiabatic argu- 
ments, one assumes that at each time the star is in a state (4.18), and the 
Hawking radiation is the result of a transition between a level |\E'(?T,tot) > 
(see (4.23)) to one of the closest levels, {"^{utot — 1) > [3]- Using energy 
conservation and (4.23), the radiation frequency satisfies 

hurad = AM{ntot, ritot - 1) ~ -j^ . (5.13) 

On the other hand the temperature is proportional to the radiating energy 
(frequency), so we have 

Th ~ huJrad = AM ~ ^ (5.14) 

in agreement with (5.12), and with Hawking results. 

In this dynamical situation, one can calculate the lifetime of the level 
l^(^tot) > [3]. This should be finite, because there is an interaction with the 
vacuum state of the radiation fields. Now this is not the Kruskal vacuum^^ 
(like in the eternal black hole case). The vacuum state is now the Unruh 
^''Known as the Hartle-Hawking [16], or Israel [15] vacuum. 
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vacuum [17]. This lifetime can be estimated to be proportional to the inverse 
of the imaginary part of the effective action (in the Unruh vacuum) , and one 
get the mass rate [3] 

^ ^ ^ ^ M (5.15) 
dt AM M2 ^ ^ 

in agreement with Hawking results, which assume a black body radiation 

rate. 

If we "extrapolate'' our results to the quantum region (n ~ 1), we can 
say that there should be a "quantum remnant" of mass M^em ~ Mp at the 
end of the Hawking evaporation [18]. But this is pure speculation because 
we ignore back-reaction as well as strong quantum effects in our model. 



6 Concluding Remarks 

In this work we used the canonical quantization approach of spherically sym- 
metric dust matter universes, first given by Lund, and apply it to the case 
of collapsing stars and black holes. The quantum states describing those 
universes are bound states and one gets a discrete spectrum. 

First let us consider some of the physical consequences of the quantized 
spectrum. One may question the collapse process itself, because if the col- 
lapsing star must satisfy (2.11) for example, then in the space of all masses 
and initial radii, only a set of measure zero satisfy it, so maybe most of the 
stars will not collapse at all? This is not the case because though the mass is 
a constant (by energy conservation), Rq (or in the general case, all the other 
geometrical quantities) can fluctuate, and one must calculate ARq/Rq. If this 
is a very small number, then the collapse is possible in a general situation. 
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Using (2.11) we have 

(6 1) 

which for astronomical objects is a very small number. For example in the 
case of our sun, we have ARq/Rq ~ 10"^°'^. This means that for astronomi- 
cal objects, the mass quantization cannot affect the classical collapse process. 
On the other hand, if we "extrapolate" our results to Plankc size objects, 
then the collapse itself may be affected by the quantization conditions. This 
may be another reason to consider stable Planck size objects (black holes)? 

The effect of the mass quantization on the Hawking radiation spectrum, 
will be mainly on very large wavelengths, A > Mbh. The black hole cannot 
radiate or absorb radiation with A > M^/i because it correspond to AM 
smaller than the distance between nearest levels. For astronomical objects 
this effect will be hard to detect. 

So we see that for "classical objects" (astronomical stars), for which n » 
1, the correspondence principle works, and the quantum effects are negligible. 

It is also quite easy to recover the classical laws of black hole thermody- 
namics: using (4.23) and A — ISnM'^ we get the first law of (Schwarzschild) 
black hole thermodynamics 

SM = 'Jisn ~ MlsA (6.2) 

And because Mbh ~ E^h we have 5*5/, ^ A and T;,/, Mp/Mbh- 

The second law is just (5.7) with the fact that AS^h > for an isolated 

system, while the generalized second law is AStot > 0, where Stot — S^h + -S". 
In the case of spherically symmetric dust matter, the infinite number 

of gravitational degrees of freedom decoupled, and each shell of dust moves 

independently. It is possible to choose the coordinates and the field variables 
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such that each shell is an harmonic oscillator. This is a simple generalization 
of the homogeneous Oppenheimer- Snyder model. 

In the case of black holes, the discrete spectrum gives the Bekenstein 
area quantization: the area of the black hole is an integer number times the 
Planck area. 

It is very easy to calculate the degeneracy of this system (of independent 
oscillators), and from it to get the entropy of the black hole. The results 
agree with the known Bekenstein-Hawking entropy: the entropy is propor- 
tional to the surface area of the black hole. Then one can use the standard 
thermodynamic definitions to get all the other thermodynamic quantities 
(e.g. Hawking temperature). 

It seems surprising that our simple model (of spherically symmetric dust 
matter) gives "enough" degeneracy, and the correct Bekenstein-Hawking en- 
tropy. One might think that in the general case (no symmetry, and a general 
matter field) the degeneracy will be much bigger, and so also the entropy, 
which would contradict known results. But this is not necessarily the case. 
One should remember that we have been able to quantize the dust system, 
because we could fix the gauge completely, which means that we correctly 
choose the coordinate system and solved the constraints. A consistent quan- 
tization of the general case (if it exists) may be achieved by a "free field 
representation" , which will be a generalization of our independent harmonic 
oscillators. If this is the case, then the general degeneracy will be quite simi- 
lar to what we have in the dust will the black hole thermodynamics. 
Maybe this is one thing that we can learn from our simple modeQ. It should 
be interesting to study some extensions of our model, and to see if our results 
will survive. Another thing that may be interesting to study in more detail, 
^^Though infinite-dimensional. 
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is the interaction between the geometrical degrees of freedom that we use 
and other fields. 
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